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Abstract 

We discuss the statistical mechanics of a gas of gauged vortices in the canonical formal¬ 
ism. At critical self-coupling, and for low temperatures, it has been argued that the 
configuration space for vortex dynamics in each topological class of the abelian Higgs 
model approximately truncates to a hnite-dimensional moduli space with a Kahler 
structure. For the case where the vortices live on a 2-sphere, we explain how local¬ 
isation formulas on the moduli spaces can be used to compute exactly the partition 
function of the vortex gas interacting with a background potential. The coefficients 
of this analytic function provide geometrical data about the Kahler structures, the 
simplest of which being their symplectic volume (computed previously by Manton us¬ 
ing an alternative argument). We use the partition function to deduce simple results 
on the thermodynamics of the vortex system; in particular, the average height on the 
sphere is computed and provides an interesting effective picture of the ground state. 

MSC (2000): 53C80, 37K65; PACS (2003): 11.27.-bd, 74.25.Bt 


1 Introduction 

One of the most challenging aspects in the study of topological solitons in gauge field 
theories is to understand their interactions, even at the classical level. At critical self¬ 
coupling, where the solitons exert no net static forces among themselves, one can typically 
describe the dynamical interactions at low speed in terms of geodesic flow for certain 
metrics on the moduli spaces of stable field configurations [1]. Exact results about these 
metrics can be obtained in some instances, and they have provided detailed information 
about the classical dynamics of solitons in this regime [2]. There is now considerable 
evidence on the beautiful geometrical fact that these moduli spaces encode a whole range 
of physical information about the underlying field theories, which goes well beyond the 
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problem of approximating the slow dynamics that brought them hrst into mathematical 
physics. 

An example illustrating how physical information can be extracted from the geom¬ 
etry of the moduli spaces is provided by the study of the statistical mechanics of a gas of 
vortices in the abelian Higgs model. Manton obtained the partition function in the crit¬ 
ically coupled (noninteracting) regime from the volumes of the moduli spaces [3]. These 
volumes can be calculated once the Kahler classes of the moduli spaces are known. At 
first, this calculation was performed for vortices on a sphere, but a similar argument can 
be used to compute the partition function for vortices living on any compact Riemann 
surface [4]. 

In this paper, the abelian Higgs model is modified by adding to the lagrangian (at 
critical coupling) an external potential that will probe the interactions of the vortices 
themselves. For weak potentials, we expect the effect of this coupling to be well described 
by the addition of a potential to the moduli space dynamics. In this setting, we calculate 
exactly the partition function for the vortex gas in a background field. Clearly, the problem 
that we consider is still less ambitious than the more physically interesting (but also more 
difficult) statistical mechanics of vortices with Ginzburg-Landau self-interaction, but it 
does provide a nontrivial extension of the study by Manton. The vortices will be allowed 
to live on a sphere with a particular axis singled out, and the potential we shall focus on is 
natural given this geometry. To obtain the statistical mechanics of the system, we will be 
making use of a localisation formula for a circle action on the moduli space of N vortices. 
This turns out to be an alternative route to obtain Manton’s results (which we recover as 
we switch off the interaction), in particular his formula for the volume of the moduli space 
of N vortices on a sphere [3]. 

2 Gauged vortices and their moduli spaces 

Given a Riemann surface S, a gauged vortex is a pair (cIa, 4>) consisting of a unitary 
connection on a hermitian line bundle C ^ T, and a section of this bundle, satisfying the 
Bogomol’nyi (or vortex) equations 


Ba^ = 0 , ( 1 ) 

HA = ^*(l-(<h,4>)). (2) 

Here, (•, •) is the hermitian structure on C, and we have fixed a metric ds'^ on S with Hodge 
star *; BA is dehned from dA = d — iA =: Da + Ba as usual through the decomposition 
provided by the complex structure of S, and Ba = dA is the curvature of dA- These 
equations are invariant under gauge transformations 

(dA,^) ^ (dA+dA,e'^<h), A G C'“(S;M) ^ aut(£, (•, •))• 

The choice of hermitian structure does not play an important role in our discussion, so 
once we hx local trivialisations for £ ^ S we use them to pull back the standard hermitian 
structure of C. 
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To a configuration {dA, ‘f’) we associate the vortex number 



Ba. 


For the cases of interest, where S is compact or is effectively compactified by imposing 
suitable boundary conditions, G Z and it corresponds to the degree of C, a topological 
invariant. Since (1) states that $ should be a holomorphic section, N is also the number 
of zeroes of <h, all having positive multiplicity. In fact, solutions of (1)~(2) are completely 
characterised by the zeroes of <I>, which can be any set of N points on T, (counted with 
multiplicity) (cf. [5]), provided that 


AttN < Vol(S) (3) 

holds [6, 7]. So the space of solutions to the vortex equations modulo gauge equivalence 
has the structure 

IT -Mn, 

tvgn 

where each moduli space of A^-vortices Ad tv is the A^th symmetric power of Z, the smooth 
2A^-manifold 

Mn = Sym^(S) := S^/6 tv. 

Complex coordinates on this (complex) manifold are usually referred to as moduli. If z 
is a local coordinate on an open set f7 C S, then the natural coordinates (zi,..., zn) on 
the cartesian product , denoting configurations with zeroes of at each z = Zr, form 
an appropriate set of moduli on Sy'm^{U) — A, where A C Ad tv is the locus on which at 
least two zeroes become coincident (and thus the STv-action fails to be free). The Zr are 
interpreted as positions of N individual vortex cores. 

The Bogomol’nyi equations above appeared for the first time in the study of the 
abelian Higgs model, a field theory for vortex dynamics defined by the functional 

Anx2{DA,^) = -\\\FA\\h + \\\DAnh - - M\h (4) 

depending on a self-coupling G M"*". The fields A and <I> here depend on a time parameter 
t G M, and the norms are taken with respect to the metric dt^ — ds^ on M x S and the 
hermitian structure on pr^T; notice Da = (^ — iAt)dt + dA(t) here has a time component, 
and we denote its curvature by Fa = d{Atdt + A) =: Ea /\ dt + Ba- It was observed [8] 
that static configurations at critical self-coupling = 1 are minima of the energy defined 
by AHi if and only if they satisfy the first-order equations (l)--(2). The (potential) 
energy of a gauged vortex with vortex number N is then ttN. Since these equations are 
easier to study than the second-order Euler-Lagrange equations of (4), one might hope 
to understand the dynamics of the abelian Higgs model from a study of gauged vortices, 
at least in the setting where the velocities are small, so that the field configurations are 
well approximated at any instant by solutions of the Bogomol’nyi equations [1]. One way 
to describe this so-called adiabatic approximation is as follows: construct from (4) an 
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action on each TA^tv by taking the helds to be solutions <h(z)) of (1)~(2) with 

time-dependent moduli (zi(t),... ,ZN{t)), and integrate over S in the norms. It has 
been proven that the resulting mechanical system does give a good description of the 
true (infinite-dimensional) vortex dynamics, even when we shift slightly from critical self- 
coulpling [9]. Upon this process of adiabatic reduction, a potential term in the held theory 
becomes a potential function for the dynamics on AIat. 

The approximated abelian Higgs dynamics on the moduli space was discussed by 
Samols [10] following work by Strachan [11], It consists of geodesic motion with inertial 
mass TT (the static energy of one vortex) for a metric grs on This metric is Kahler 

with respect to the complex structure on A^jv induced by the one on S. We call it the 
metric on A^at, since it is obtained from the natural norms (4) on the space of 
(covariant) derivatives of pairs (Da,^)- It is described by the closed (l,l)-form 



N 

QrsdZr A dZg 

r,5=l 


i f db \ 

- ^ ( ^'^{Zr)5rs + 2^ j dZr A dZg. 

r,s=l '' ^ 


( 5 ) 


Here, Zr are the moduli associated to a complex coordinate z for which ds'^ = {z)\dz\‘^. 
The quantities br{zi ,..., zjy) are dehned as follows. One can combine (l)-(2) into a single 
equation for the gauge-invariant quantity h := log(<h,<I>), 


d^h 

dzdz 


+ n‘^{z){l-e^) 


N 

dvr 5{z — Zr), 

r=l 


( 6 ) 


where 6 is the Dirac delta-function. A solution h{z; zi,..., zn) to (6) has an expansion 
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h{z) = log \z — Zr\^ + ttr -\ - br{z — Zr) -|- H- br{z — Zr) + 




\Z - Zr\ + 


( 7 ) 


about a simple zero Zr of <I>. It is a remarkable fact that only the linear coefficients hr in 
this expansion appear in (5). It is clear that (7) only makes sense in a neighbourhood of 
the vortex Zr which does not contain any other vortex positions. In fact, (6) implies that 
we can write [10] for each r = 1,..., A 


N 

hr{zi,. . . ,Zn) = 

s^r 



+ brizi,. . . ,Zn), 


( 8 ) 


where hr are smooth on the coincidence locus A. In an expansion about the position 
of a vortex with multiplicity n < N, the logarithmic term in (7) has a prefactor n. 
The local moduli Zr cannot be used to describe a neighbourhood in A4 at of such a vortex 
configuration, and so the coefficients in an expansion equivalent to (7) cannot be expressed 
as functions of the Zr- 

Samols’s formula (5) still does not give the Kahler form explicitly, since the nontrivial 
quantities br are specified in terms of unknown solutions to (6). Extracting concrete 
information about the moduli space metrics is still a nontrivial challenge. One may feel 
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that work in this direction must inevitably rely on a numerical study of equation (6); 
however, analytical results have been derived in particular cases using approximations of 
some sort [12, 13], integrability [11], or even a remarkable argument involving T-duality in 
string theory [14]. In the following, we shall obtain further analytical information about 
the moduli space metrics when S is a 2-sphere. 


3 Vortex dynamics in a background potential 

In this paper, we would like to extend the abelian Higgs model to include a coupling 
with a background potential, which we define to be any smooth function / : S ^ M. 
To introduce the coupling term, we first define a vorticity 2-form v{Da,^) on S by the 
equation (cf. [15]) 

{dj + Fa) a dt =: V a dt, 

where j is the gauge-invariant supercurrent 1-form on R x S 

j := lm{(j),DA(l)). 

The coupling to the potential / we propose is given by adding to the potential energy 
term of the action functional (4) the interaction term 
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fv 


(9) 


where fj? is a coupling constant. 

The Euler-Lagrange equations for the modified action 

An^2 {Da, j f v{Da, d>) A dt (10) 

are as follows. We obtain the same Gaufi’s law as for the abelian Higgs model (4), 

d * Ea = 2 *lm Dt^ dt) (11) 

(where Dt := ^ — iAt), which is to be regarded as a constraint on the space of fields 
determining At- The dynamical equations are 

d*BA- = - * Im (4>, dA^) - $) + 1) df, (12) 

ot 
\2 

= y ((^, - 1)^ + 2i/i2 * {df A dA^), (13) 

where □ := {Dt)‘^ — * dA * dA is the covariant d’Alembertian on pr^T ^ R x S. As 
expected, there are new terms (proportional to /A) adding to the current in Ampere’s law 
(12), and also to the potential in the nonlinear Klein-Gordon equation (13). 

From our discussion in section 2, we know that solutions to the vortex equations (1)- 
(2) solve the equations of motion (11)-(13) above in the static case (and setting At = 0), 
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provided = 1 and = 0. When the couplings are perturbed slightly away from these 
critical values, we expect a slow-moving solution to have a best vortex approximation, and 
that we can follow its evolution under the dynamics defined above. A detailed analysis 
of the case ~ 1, = 0 has been carried out by Stuart in [9], where distances to best 

vortex approximations were estimated (in terms of A^ — 1 and initial errors) after finite¬ 
time evolution, and shown to be controlled for evolution times of order O (|A^ — 

In this paper, we make the natural assumption that an analogous result holds for A^ = 1 
and ~ 0. This assumption can be regarded as a motivation to the study of the dynamics 
on the moduli space of gauged vortices with a certain (and rather natural) background 
potential. 

The coupling (9) satisfies a number of desirable properties. One of them is that, if 
we take / to be a constant, it reduces to a constant potential on the moduli space. To see 
this, we use (1) and (2) to rewrite [15] 

V = dA<^) + (1 - (d>, <^))Ba 

= -i{dA^, dA^) + 2 {*Ba)Ba, 

which is precisely twice the energy density for solutions of the Bogomol’nyi equations at 
critical self-coupling. Thus 



Notice that df = 0 in this situation, so the terms proportional to fj? in the field equations 
(12) and (13) vanish. This is a rather degenerate case of our model (10), but the fact that 
the adiabatic approximation leads to sensible results at this level is already reassuring. 

For the rest of the paper, we shall restrict ourselves to vortices living on a 2-sphere 
of radius i?, S = In this context, we shall illustrate another natural property of the 
interaction term (9) in section 4.1. 


4 Localisation and the partition function 


The exact results we want to derive refer to the case S = 5^, on which z shall denote a 
stereographic coordinate. The metric on has Kahler (volume) form 

2iR^ 

cosi = 


and the constraint (3) reads 

The moduli spaces in this case are 


« (l + kP)2 

i?2 > N. 


dz A dz 


Mn = Sym^(5|) ^ CP^, 
equipped with the Kahler structures 

N / 


to = I 


‘E 


R^Srs dbg , 

(J Zr- 


V(1 + krp) 


(14) 


(15) 


(16) 
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One way to understand (15) is as follows. Let C AItv be the locus where precisely 
j vortices are at z = oo. Then Vj is parametrised by the (unordered) positions of the 
remaining N — j vortices on C, which are unambiguously specified by the coefficients of a 
monic polynomial of degree N — j having these positions as roots. Hence Vj = . The 

way in which these 2{N — j)-cells are glued together in AItv is determined by attaching 
maps corresponding to letting one vortex go to z = oo at a time. This yields precisely the 
description of CP^ as a CW-complex. 


4.1 Rotational symmetry 

For most of our discussion, we shall restrict our attention to the potential 


(17) 


This potential is very natural if we assume that there is a special axis on the sphere. In 
fact, (17) is the simplest nontrivial circularly symmetric function on 5^, in the sense that 
other potentials with circular symmetry can be expanded as power series in /. 

It is easy to check that (17) is a hamiltonian for the circle action of rotations around 
the axis of associated to the stereographic coordinate z, 




There is an induced circle action on AIat 


Sym'^(5^) with generator 



(18) 


which extends to a smooth vector field on AIat. Rotational symmetry on implies [16] 

N 

(^Zrbr — Zrhr) = 0, (19) 

r=l 


and one can use this equality to show that the Kahler structure (16) is preserved by the 
one-parameter group generated by (18). Since cu is closed, it follows from 


0 = £^uj = i^{duj) + = d{i^uj) 

and H^{CF^) = 0 that there is also a hamiltonian for the circle action on AIat, which can 
be computed to be 

^ / 1 _ U |2 \ 

J = -Mr + l)j . 

r=l '' ' 2 

Notice that this formula assumes that all the vortices are separated, but it does extend to 
a smooth function on the whole of AItv- 
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Hamiltonians of circle actions are defined up to a constant. However, this constant is 
fixed if the circle involved is a one-parameter subgroup of a Lie group with discrete centre 
and hamiltonian action on the symplectic manifold, and we demand that the hamiltonian 
should be a component of the corresponding moment map. In our case, the circle action 
extends to a hamiltonian action of Iso(5'^) = SO(3) on Mn- The choice of constants in 
(20) can be checked [16] to be consistent with the moment map Mn 5o(3)*. The same 
can be said of / in (17) in relation to the symplectic structure u)g 2 in (14). 

A natural question to ask at this point is whether our coupling behaves well with 
respect to moment maps. We can regard the adiabatic reduction of couplings of the form 
(9) as defining a linear map between Lie algebras 


: 


f 


fs%f^ 


( 21 ) 


for each N (where the Lie bracket on each side is the Poisson bracket defined by the 
symplectic structure). One may hope that this map preserves the Lie algebra structures, 
and that it relates corresponding components of the SO(3)-moment maps on each of 
and Mn- The following proposition shows that this is indeed the case. 

Proposition 4.1. For each N < B?', there is a commutative diagram of Lie algebras 


C^{Sl) 


so(3) 


n 


where the diagonal arrows denote the dual maps to the moment maps ^ so (3)* and 
Adw ^ 50(3)*, respectively. 

Proof. The existence of the moment maps is guaranteed by the vanishing of the Lie algebra 
cohomology group L7^(6o(3); M), a consequence of simplicity [17]. 

Given the linearity of TZ in (21), the proposition will follow if we check the commu¬ 
tativity of the diagram on the generators of so (3). 

We start by sketching how to obtain 

J = 7^(/), (22) 

where / is given by (17). This calculation is paradigmatic of the process of reduction to 
the moduli space. Suppose that the vortices are all separated, and work first on the subset 

r n ^ 

C,:=lz£C:\z\<-\-[j B,{zr) C 3% 
where e is taken small enough. Then (6) can be used to write on 
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and 


fv = 


I i£(i - ( S I + i£(i + 


d 


d fdhV d 


2R^^\dz' 


dz \dz J 


dz' 


dz 


A dz. 


Using Stokes’ theorem and the estimates 
d^h dh 


dzdz 


■(1+|Z.|2)2 


+ o(e), 


dz e 2 


as e 


for z G dBf:{zr), which follow from (7), we do obtain (22) after taking e ^ 0. 

To proceed, we can either repeat the calculation for the other generators (as in [16]), 
keeping track of the constants to preserve the so (3) algebra, or change integration variables 
and apply the appropriate rotations to the two sides of (22). The second procedure 
becomes straightforward once we observe that the hr transform as [4] 


{T*hr){zi, ...zn) = 


1 .. ^ T"{zr) 


under any holomorphic T G Iso(S). This equation is readily obtained from the expansion 

( 7 ). 

□ 


It should be noted that TZ does not preserve the structures of C°°{Sj^, u)g 2 ^) and 
C°°(J\Ajv,uj) as Poisson algebras. 

Since the moduli spaces are compact in our case, the circle actions we are interested 
in must have fixed points (the zeroes of ^). It is easy to see that these are exactly the 
+ 1 points pj G MIV (with j = 0,1,... ,N) describing configurations of j vortices at 
z = 0 and N — j vortices at z = oo, which will be fixed by a rotation of all the vortices 
around the axis through 0 and oo. We remark that J is a Morse function on AIat, with 
critical set 

Crit(J) = {po,Pi,... ,Pn]- 

In the next section, it will be useful to understand the circle action in the neigh¬ 
bourhood of the fixed points. On each tangent space Tp.Aijy = C^, the action linearises 
to a complex Wdimensional representation of the circle group; this in turn decomposes 
into N 1-dimensional representations of U(l), and each of them is uniquely specified by 
its weight k G Z: 

: C ^ tGR, C e C. 

Thus to each of the fixed points pj we can associate N weights kj^£, 1 < i < N, well 
defined up to order. The product of the weights at each fixed point 

N 

e(Pi) := n (23) 

e=i 

is a local invariant of the circle action; it does not depend on the choice of coordinates on 
the moduli space. We shall make use of the following result; 
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Lemma 4.2. For the circle action generated by (20) on 


e{Pj) = - j)\, 0<j<N. 


Proof. In terms of the moduli Zr, the circle action is simply given by 




^27rit, 


t e 


( 24 ) 


However, all the pj except po lie in A, where the coordinate system dehned by the Zr 
becomes singular, and so we must introduce other coordinates to compute the weights. 
We hx j and arrange the vortex labels such that vortices 1,..., j are at z = 0, and vortices 
j + 1,N are at z = oo. (We are allowed to do this since the vortices in each cluster 
are to be thought of as interchangeable, but vortices belonging to different clusters have 
separate identities.) Now we introduce 

Uj^r ■= Sr 1 < r < j, 

Vj,r ■= (z“,_\,..., z)^^), I <r <N -j, 

[?1 

where s). denotes the rth elementary symmetric polynomial in j variables. 

Clearly, {uj^i ,..., Ujj,Vj^i ,..., Vj^N-j} is a centred local coordinate system at pj G M.n, 
and we can also use it as a coordinate system on Tp.Mjsf. From (24), we find that the 
linearisation of the circle action is described by 


and therefore we obtain in (4.2) 

e{pj)= (^n 


^27Tirt^ 






/N-j \ 

\t=l / 


□ 


4.2 The partition function 

The dynamics on the moduli space of N vortices is defined by the lagrangian 

N 

L = - ^ Qr^Zi, . . .,ZN)Zrts - M^J(zi, . . . ,ZAr), 
r,s=l 
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where grs are the coefficients of the metric on and vr is the mass of a single vortex. 
In the canonical picture, we describe the dynamics as a hamiltonian system on the 4iV- 
dimensional manifold T*M]sf equipped with its canonical symplectic structure 


^can 


1 

2 


N 

'^^{dZr A dWr + dZr A dWr)- 

r=l 


Here, the Wr are complex coordinates for the fibres of 
moduli Zr'. 


Wr 


dL 

dZr 


N 

^ ^ 9rsZs- 
s=l 


Mn, conjugate to the 


The dynamics is generated by the hamiltonian 



N 

. . .,ZN)WrWs + g^J{zi, . . .,Zn), 

r,5=l 


where denote the entries of the inverse to the matrix of coefficients of the metric. 

According to the canonical formalism for classical statistical mechanics, the partition 
function of the vortex moduli space dynamics is given by the Gibbs formula 


Z = 


OJ. 


2N 


{2'Kh) 


2N 


It*Mn 


{2Ny: 


(25) 


The prefactor to the integral is Planck’s constant 27rfi raised to the power ^ dimR T*M-n = 
2N, and T is the temperature. We are normalising Boltzmann’s constant to unity. 

In parallel with the calculation in [3], we find that the integral (25) factorises as a 
product of a gaussian integral along the fibres, which can be readily calculated, and an 
integral over the moduli space with the Liouville measure associated to (16): 


Z = 



N\' 


(26) 


The integral remaining still looks very complicated, but we shall show that it can also be 
computed exactly, using localisation in symplectic geometry. The main tool we will use is 
the following version of a famous result by Duistermaat and Heckman [18]: 


Theorem 4.3 (Duistermaat—Heckman formula). Let {M,uj) be a (2n)-dimensional 
compact symplectic manifold with a hamiltonian circle action generated by a Morse func¬ 
tion iL : M ^ M. Then 



E 

p&Cnt{K) 


^tK{p) 
T^e{p) ’ 


(27) 


where e{p) denotes the product of the weights of the linearised action at the critical point 
p, and r is a formal parameter. 
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A streamlined proof of this theorem can be found in section 5.6 of reference [19]. 

In our context, taking M = M.n, K = J and r = and making use of 

Lemma 4.2, the equality (27) takes the form 




N 


iMp 


OJ 

ivT 


N 

j=0 




(28) 


Hence, to compute the integral in (26), we only need to evaluate the potential J at the 
critical points pj £ Mn- 

We shall make use of the spherical symmetry of the problem to determine the con¬ 
tribution of the br terms in the formula (20) to J{pj)- Suppose that the Higgs field 4> has 
a zero at z = y of order j and a zero at z = — 4 of order N — j, a vortex configuration 
corresponding to a point of AIat that we shall denote by Pj{y)- Then h = log(<h, 4>) must 
satisfy 

d‘^h i?^(e'‘-l) _ ( 29 ) 


= jirdiz — y) + (N — j)Ti5 \ z — - 

y 


dzdz (1 -|- \z\‘^)‘^ 

This equation leads to the following expansions for h: 

1 


1 - 


hiz,y) =jlog\z - y\ +a+{y) + -b+{y){z - y) + -b+{y){z - y) H- 


around z = y, and 

h{z,y) = {N -j)\og 


1 

z + - 

y 


+ «-(y) + o&-(y) 2: + - +7;b-{y)[z + - + 


around z = —1. We need to calculate (y). Extending an argument in [3], we explore the 
fact that h(z, y) must be a function of the spherical (and thus also of the chordal) distance 
between the points of coordinates z and y on S^. The square of the chordal distance is 
given by 

4i?2|^-y|2 


(l + |z|2)(l + |y|2)’ 


and so h can be expanded as 


h(z, y) = j log 


\z - y|" 


(l + |z|2)(l + |y|2) 


+ c+ -|- dj^ 


jz - y'\‘ 


(l + |z|2)(l + |y|2) 


+ 


(30) 


For y = 0, this yields 


h(z,0) =y logjzj^ + c+ -k {d+ - j)\z\^ + 


substituting this into (29) with y = 0 and small z, and looking at the zeroth order term 
in jzp, we conclude that 

d+=3-R\ 
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Now taking y arbitrary, we rearrange (30) as 


Hz,y) = j^og\z-y\ + c+- 2jlog(l + |y| )- 

.,2 


2^ jy f X jy i- -N 


+ 


jy"^ ( x2 , jy (- -n2 

\Z y) “1“ L.i9\9 y) 


l + |y|2^ 1 + |2/|2' 


(1 + | 2 /P) 2 ' 


to read off 


(i + |yP)2 


^+(y) = , n.|2 - 


(i + |yP)2 


l^-2/l + 


We proceed similarly for b-{y): writing 
h{z,y) = (A^-j)log- 



Z + } 

2 


y 





1 + 


+ c_ + d— - 


( 31 ) 



Z + }- 

2 


y 



(i + H^) (' + #)' 


+ 


(32) 


we can calculate 

C-=N -j-R^ 

by taking y = oo, substituting in (29) and looking at the zeroth term in the large \z'\^ 
expansion; then rearrange (32) as 


h{z,y) 


{N -j) log 


1 

z + - 
y 


2 


+ 

+ 


{N-j)y I 
l + |y|2 ' 

(N -j)f 
(i + |yP)2 



- 2(A^ - j) log 
, {N-j)y f 

i + |yP V 


1 + ll/ p 
|y| 

1 

z + - 

y 


+ C_ 


+ (l + |y|2)2 y) 


{i + \y?? 


1 

^ + - 
y 


2 

+ ... 


and read off 


b-{y) 


‘^{N -j)y 

l + W 


(33) 


To proceed, we must interpret the formula (20) carefully to deal with vortex clusters 
such as the configurations Pj{y) G Mn- When any number of vortices become coincident, 
the smooth part br of their individual br coefficients tends to the linear coefficient in the 
cluster expansion of h about the position of the cluster, but the singularities in (8) must 
be treated with some care. Although these singular parts diverge separately, they always 
yield a finite contribution in a Sjy-invariant linear combination over the individual vortices 
— for example, they cancel mutually in a term like ^r- In (20), these singular parts 

are multiplied by the vortex positions, and so in the clustering (zi,..., zn) Pj{y) they 
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give a contribution 
N 


3 


N 


-2'K'^Zr ^ 

r=\ s^r 

same cluster 


2:r — 2:c 


2:^ — 


E E 


N 


-4^ EE 

r=l s^r 

-^'(ee+ee 

\ r=l s<r r=j+ls<r 
/ j N 

-4- EE+ E Eb 

\ r=l s<r r=j+l s<r, 


r=j+l s^r 
Zr 


Zr — Z. 


Zr — Zg Zg — Zr 


= —47r 


j(j-l) I {N-j)iN-j-l) 


2 2 
= -27r{iN - jf + f - N) 

to J{pj{y))- The remaining part of J{pj{y)) is 


(34) 


N 


2vrE 


r=l 


,1 — \Zr 

\ -\-\Zr 


— (Zrbr + 1) = 27ri?^ j 


2 A-\y\' 


1 + I2/F 


+ iN-j) 


1 - 




1 


-27r \^j{yb+{y) + 1) + {N - j) ) b-iy) + 1 

= 27r(2j-iV)(i?2 _ AT)^ “ 


'I + I2/P 

+27r{{N-jf+f-N), 

where we used (31) and (33). Adding (34) and (35), we obtain 

i-|yp 


(35) 


and hence 


iMr 


J{pj{y)) = 

27r{R^ 

= J{Pj{0)) = 

2tt{R^ 

(28), we find 

, A 

N 

m ■ 

= E( 


1 + |2/|^ 


j=0 


rpN ^-2-Kf?{R^-N){2j-N)/T 


(36) 


(37) 


This equation can be interpreted as an equality in the formal power series ring 
which is equivalent to an infinite number of identities over the reals. The hrst N nontrivial 
identities are 


N 




(38) 


and they must be true for consistency. But they are implied by the following technical 
lemma. 
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Lemma 4.4. For all N G'N 


{-ly fN\ (N 


m \j)\2 V 


0 if Q<k<N -I, 
1 if k = N. 


Proof. We start from the equality 




obtained by successively acting on the binomial expansion with the Euler operator x-^. 
Setting X = 1 yields 


We use this to write 




N\ fN 




^ ^ / at\ /i \ / Ar\ ^^ 


j=o e=o 


E(-i)' 


k\ fN 


k-t N 


E(-F ” / 


E(-F 


k\ fN 




X = 1 

It is clear that (39) is zero whenever 0 < k < N, since the differential operator ( y ~ 3;^)^ 
will not annihilate enough (1 — x) factors before x ^ 1. However, ii k = N one obtains 




= (1-* 


= N\ 


making use of = 1. 

Notice that Lemma 4.4 yields yet another identity from (37); 


Vo1(M7v) := 


(47r)^(i?2 - N)^ 

W\ 

(47r)^(i?2 - N)^ 

W\ ■ 


(-1)^' (N 


m-jv- V2 


15 



This is precisely the formula found by Manton for the volume of the vortex moduli space 
in [3], using a more direct argument involving the cohomology of CP^. Equation (40) 
provides a nontrivial check of our calculations. 

Beyond the identities (38) and the formula (40) for Vol(A4Ar), our localisation argu¬ 
ment yields an infinite number of integrals over the moduli space: for m G N 



(—1)^ 


{2Tr{R^ 


N){2j-N)f^’ 


(41) 


Notice that both sides of this equation vanish when m is odd: The left-hand side from 
reflection symmetry on and Proposition 4.1, and the right-hand side from 


N 


E 


(-ly 

JKN-jV. 


\ N-\-2n—l 


Vn G N, 


which follows from substituting j hy N — j in the sum. For m even, (41) yields new 
quantitative information about the metric on Mn- We note in passing that even the 
vanishing integrals (m odd) have interesting content. For example, using the result for 
m = 1, we find from (20) that 

/ ^2:r6,.(2:i,...,2:Ar)— = -Vol(A4Ar); 

JMn f—i ■ 

a similar argument for rotations around z = 1 and z = i yields 


f ^ ui^ 

/ y'6,.(2:i,...,2:jv) —= 0 


which in turn also leads to 



UJ 


N 




by a property of the functions br analogous to (19) [16]: 


N 

^ ^^ (2zr z^br 6j.) = 0. 

r=l 


Using the result (41), together with (40), we can obtain the integral over AIat of 
any power series in J. Such power series, when convergent, give analytic functions on AIat 
which are invariant under the circle action generated by (18). 


16 



Finally, we can use (37) to compute the partition function (26) of the vortex gas in 
the backgound field (17) to be 


Z = 


( V p2ni,-^N{R^-N)/T ^ (~^)^ -4nu^(R^-N)i/T 


J_ f ^ ^ 2■K^?N{R^-N)/T (1 _ -i^T^?{R^-N)/T\ 

m ) V ) 

2TT^i‘^{R^ -N) 


N 


N\ \hiiJ 


(42) 


5 Thermodynamics of the vortex gas 

The Helmholtz free energy of the vortex system can be computed from (42) as 


F = —T log Z ~ —NT log sinh 


IJ?{A — AttN) 


2T 


- log V + 2 log 

% J 


where we made use of Stirling’s approximation IniV! ~ Nh\N — N, and we introduced 
the area of 5^, A = The entropy is given by 


dF ( 

S = = N ( logsmh 


^^{A — At^N) 

2r 


+ log 




n'^iA — AttN) p^iA — A-kN) \ 

--coth-—- 

2T 2T J 


Both these quantities turn out to be nonextensive, due to a nonlinear effect produced by 
the interaction with the external potential. 

The interaction is controlled by the coupling fi^; at small coupling, keeping A and 
N finite, we can approximate the hyperbolic functions to first order as sinhx — X smd 
coth x — ^, leading to the same results found by Manton in the absence of interaction [3]. 
In this noninteracting setting, both F and S become extensive in the thermodynamical 
limit of yl ^ oo, N oo and constant density n = ^. The pressure P = — ^ of the 
system in this regime can be readily computed and yields the equation of state [3] 


P{A - AttN) = NT. 


(43) 


This is a particular limit of the van der Waals equation, known as a Clausius equation 
of state. It holds more generally on any compact Riemann surface [4]. The fact that the 
factor A —AttN appears in (43) can be interpreted as an interaction among the vortices [3]: 
each vortex effectively occupies an area of Att (consistently with (3)), hence N coexisting 
vortices have an area available for their motion which is a reduction of the area A of the 
sphere by V x 47r. The virial coefficients associated to (43) are found to be all constant 
and equal to powers of dvr: 

OO 

PA = iVr^(47r)V; (44) 

^=o 
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this virial expansion is reminiscent of the one for a gas of hard particles of finite size in a 
one-dimensional box [2]. Thus one might be tempted to think of the vortices in effective 
terms as rigid discs of area dvr moving on the surface. But this picture already fails at 
first order in the expansion (44); the gas of hard disks would have first virial coefficient 
Stt [20], which is twice the coefficient of n in the power series in (44). In fact, a crucial 
difference between the gas of vortices and a gas of hard discs is that the shapes of the 
regions where the vortex density is concentrated become very different in the two cases 
whenever two or more particles come close together, and we expect the equation of state 
to be extremely sensitive to this. 

Now we shall consider the statistical mechanics of the interacting regime ^ 0). 
From our partition function (42), we can compute the thermodynamic average value of 
the observable J G C°°{M.n) given by ( 20 ): 


{J) 


1 

Z 


iT*Mn 

-NT 


, ,2N 

j ^-H/T ^can 

{2N)l 


A — AnN , u?{A — AttN) 
coth- ^ ^ 


2T 


2T 



(45) 


Recall that J is related to the height function on the sphere X 3 := G 

by J = TZ[Rx^) (cf. (22)). Thus we propose to interpret the quantity J/{2ttNR) as an 
observable giving the height on 5^ of configurations of N vortices, and we shall denote it 
by X 3 . We find 


(^ 3 ) 


jJ) 

2'n:NR 



^coth y 



(46) 


where 

Ijl^{A — At:N) 
X ■- ^ 


(47) 


The dependence of the average height on the parameter y is plotted in Figure 1. 

To interpret the meaning of (46), consider a simple model where the vortex density 
is supported and homogeneously distributed on a spherical disc of area AttN and centred at 
the minimum z = 00 of the potential (17). The height of the points at the boundary 
of this disc is computed as 


27ri? 


dx3 = AttN x\ 


'-R 




where we made use of Archimedes’ hat-box theorem. In this model, the average height for 
the vortex density is then given by 


^^3)model _ 


^„max 

X3dX3 


f-R dX3 



(48) 
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Figure 1: The thermodynamic average (xs) of the height observable in the gas of N vortices as a 
function of the parameter x- 


Thus we find that 

hm (X 3 ) = 

X^oo 

In other words, the effective shape of an Al-vortex distorts to a spherical disc (of area 4ttN 
and centred at the minimum of the potential 2 : = 00 ) as x becomes large. This limit is 
typically achieved if the area A — inN available for the dynamics is large, the coupling 
is large, or the absolute temperature T is low — cf. (47). Notice that in the picture 
where each vortex is approximated by a rigid disc of area dvr, the average height of the 
density distribution would be higher than (48). The picture of the N-vortex as a disc of 
constant density (localised at the bottom of the potential) provides a good description of 
the ground state of the system. 

To describe a circularly symmetric distribution of vortices on the sphere, it is natural 
to introduce a vortex density function p : [—R, i?] ^ M, defined on the interval of heights 
on and normalised as J^j^pdxs = 1. This function has parameters A, T and N, 
and it could in principle be reconstructed as a Fourier-Legendre series from the partition 
function (42) if we could compute the reduction to the moduli space of all the powers of 
the height, 77(/™) for m G N. We can write 

(J) 

:= {i3) = 2 ^. 

which gives just the first Fourier-Legendre coefficient. The difference in pressure AP 
between the highest and the lowest points of the gas on the sphere depends only on the 
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trivial zeroth order coefficient: using the relation 

VP + NpVf = 0 

for the pressure of a fluid of N particles subject to a potential / at equilibrium 
analogous to the problem of fluid motion in a constant gravitational field, cf. [21] [ 
find 

/ Ft pFt pi 

VPdx3 = - Npix^) — (^?RxP dx3 = -p^NR < 0. 

-R J-R 9X3 

6 Discussion 

We have been able to calculate the partition function of a gas of critically coupled abelian 
Higgs vortices interacting with an axially symmetric background potential on a sphere. 
When we switch off the interaction, we recover Manton’s partition function that gives 
physical insight on the Liouville volume of vortex moduli spaces Mn- Our study yields, in 
addition to Manton’s formula (40) for Vol(A4Ar), an infinite number of nontrivial integrals 
(41) over Aijy- These are additional data about the geometry of the moduli spaces, and 
they are all encapsulated by our partition function. As an application, we computed 
the thermodynamic average height of the gas of N vortices, and the result we found is 
consistent with the effective picture of the ground state as an Wvortex localised at the 
bottom of the potential as a spherical disc of constant density and area AttN. 

Our analysis was essentially an application of the Duistermaat-Heckman localisation 
formula for a natural circle action on in which the symplectic structure of the 

moduli space is a crucial ingredient. There is an alternative model [22] for dynamics of 
Ginzburg-Landau vortices with a Schrodinger-Chern-Simons kinetic term, for which Ai n 
(not T*Aijsr) plays the role of phase space in the adiabatic approximation; in this context, 
the Kahler form lo appears naturally as a symplectic structure [16]. Our work illustrates 
that the symplectic point of view can also be fruitful in the study of the abelian Higgs 
model. 

We have already noted that our formula (41) can be applied to calculate integrals 
of general circularly symmetric functions on the moduli space. One may therefore hope to 
study the interaction of the vortices with any symmetric potential / on — or perhaps 
even an SO(3)-invariant intervortex interaction modelling the Ginzburg-Landau potential 
at / 1 to some degree of approximation, which would be an obviously interesting 
extension of our work [23]. Analytical results about the Ginzburg-Landau potential (for 
arbitrary N) are already available [24, 25], but they refer to the situation where the vortices 
are well separated on the plane. A treatment of the interaction to include the interesting 
effects near clustering configurations on the moduli spaces will almost certainly need to 
use some numerical input [26]. It is believed that abelian Higgs vortices even slightly away 
from critical coupling should satisfy a realistic equation of state such as the van der Waals 
equation, which accounts for phase transitions. Progress in this direction would shed light 
on the phenomenology of thin superconductors. 


(this is 
;25), we 
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